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ABSTRACT 

We searched for radio pulsations from the potential neutron star binary companions to subdwarf B stars HE 0532-4503, HE 0929- 
0424, TONS 183 and PG 1232-136. Optical spectroscopy of these subdwarfs has indicated they orbit a companion in the neutron star 
mass range. These companions are thought to play an important role in the poorly understood formation of subdwarf B stars. Using 
the Green Bank Telescope we searched down to mean flux densities as low as 0.2 mjy, but no pulsed emission was found. We discuss 
the implications for each system. 



1. Introduction 

The study of millisecond pulsars (MSPs) enables several types 
of research in astrophysic s, ranging from binary evolution (e.g. 
lEdwards & BailesL l200lh . to the potential detection of back- 
ground gravitational radi ation using a large set of pulsars with 
stable timing properties ( Jaffe & Ba cked 12003 ). Furthermore, 
high-precision timing provides opportunities for fundamental 
physics research: by measuring and modeling pulse arrival times 
from binary systems, one can constrain neutron star masses and 
equations of state (Demorestet al., 2010), or study strong-field 
general relativity (ITavlor & WeisbergLll989l) . 

Finding such millisecond pulsars is done through wide-area 
surveys or through targeted surveys, each with specific advan- 
tages: with limited telescope time, a directed search can afford 
to dwell on a specifically interesting part of the sky for longer 
than an undirected search. A directed search is therefore more 
sensitive to faint pulsars in specific locations than an undirected 
search would be that only incidentally scans over such a location. 
An undirected search however covers a larger area. Directed 
searches have tar geted sources associa ted with MSPs, such as 
globular clusters ( Ransom et all l2005h . steep sp ectrum sources 
("Backer et al., '1982'), unidentified Fermi sources (Ransom et al., 
|2011) and low mass- white dwarfs (van Leeuwen et al., 2007; 
lAgiieros et all l2009i) . Here we report the results of a directed 
search for radio pulsations from four short-orbit sub-luminous B 
dwarf (sdB) binaries ide ntified through spec troscopy as possibly 
contammg neutron stars (iGeieret aUl2008l) . 

Sub-lumi nous B dwarfs are s ome of the most abundant faint 
blue objects (iGreen et al 1 Il986h . They are thought to be light 
(about 0.5 Mq) core helium burning stars with very thin hy- 
drogen envelopes. After the helium burning in the core stops, 
sdB stars evolve directl y to the white dwarf cooling sequence 
(cf. a recent r eview by Hebeil 20091) . In their survey of sub- 
dwarf B stars iMaxted et alj (l200l'l) find that 2/3 of field sdB 
stars are in short-orbit (P <~ 10 days) binaries. Taking into ac- 
count the insensitivi ty of their survey to longer-period binaries, 
iMaxted et alj (1200 lb conclude that binary star evolution is fun- 
damental to the formation of sdB stars. 



Several such binary formation channels have been hypoth- 
esized. For an sdB to form, a light star must lose most of its 
hydrogen envelope and ignite helium in its core. In these binary 
systems, sdB stars can be formed through phases of Common 
Envelope evolution where the envelope is ejected, or through 
stable Roche Lobe overflow stripping the donor star of its hydro- 
gen envelope. Binary population synthesis models for the above 
scenarios in the case of white -dwa rf companions w ere recently 
compaied by|Hanetal] (l2003h and lHu etal ] (l2007l) . 

Specifically interesting for finding new millisecond pul- 
sars, is the channel that leads to short-period binary systems 
containing an sdB star and a neutron star or a black hole 
('Podsiadlows ki et aI.Ll2002l:lGeier et"ani2010l) . Here, the binaiy 
system contains a star that is massive enough to evolve into a 
neutron star or a black hole and a secondary, lighter star that is 
the progenitor of the sdB star The secondary is in a wide or- 
bit with a period of about 20 years. The system undergoes two 
phases of common envelope evolution, a supernova and a short 
X-ray binary phase. A first common envelope phase starts soon 
after the primary becomes a red supergiant and starts overflow- 
ing its Roche Lobe. During this first common envelope phase 
the orbit of the binary shrinks. After the subsequent supernova, 
the primary leaves a neutron star or black hole. This neutron 
star remains visible as a normal radio pulsar for on average 
10^"^ years, and turns olf. Once the secondary evolves olf the 
main sequence the second phase of mass transfer starts when the 
secondary begins to overflow its Roche Lobe. The system un- 
dergoes a s hort X-ray binary phase, in which the neutron star 
is recycled dBhattacharva & van den HeuvelL[T99Tl) . The second 
common-envelope phase starts shortly thereafter, further shrink- 
ing the orbit of the binary and dissipating the envelope of the 
secondary. The secondary, which is now mostly stripped of its 
hydrogen envelope, ign ites helium in its core and continues its 
evolution as a sdB star dPodsiadlo wski et aP. '2002^ iGeier et all 
2010). From population modeling, Pfahl et al. (2003]) find that 
about 1 % of sdB stars evolve as outlined above and orbit a neu- 
tron star or black hol e. In such syst ems the secondary is an sdB 
star for ~ 10 years (lHeberll2009h . before evolving to a white 
dwarf. 



2 



Thijs Coenen et al.: A search for radio pulsars in four subdwarf B star binaries 



While low-mass companions (white dwarfs, main se- 
quence stars) to sdB stars have been optically confirmed (e.g. 
[Oeieretal., 2010), the unambiguous identification of a high- 
mass companion in the form of a radio pulsar would provide 
further constraints on the relative likelihood of diff'erent binary 
evolution models (as was similarly done for com panions of low- 
mass white dwarfs in Ivan Leeuwen et all '2007). Furthermore, 
timing of a possible MSP companion will likely provide an sdB 
mass measurement, either through the modeling of the orbital 
parameters, or more dire ctly and precise l y thro ugh a Shapiro 
delay measurement (e.g. iFerdman et al.L |2010|) . As millisec- 
ond pulsars can turn on shortly after the cessation of accretion 
(lArchibald et al.L 12009) and can shine for more than the ~ 10'^ 
year sdB star lifetime, MSPs will be active throughout the sdB 
star phase. Thus, non-detections of radio pulsations from these 
systems could mean the absence of any neutron star, the pres- 
ence of only a very weak radio pulsar, or of a brighter one that 
is beamed away from Earth. Non-detections in a large enough 
sample of sdB stars eventually provide statistics on sdB forma- 
tion channels. 

In Section |2] we describe how sources were selected; in 
Section [3] our observing and data reduction setup is outlined. 
Section |4] contains our results. In Section |5] we discuss these 
results and compare to recent optical results on some of our 
sources. 

2. Source selection 

Our so urce selection was based on the report by iGeier et aP 
(1200 8') of the detection of the first four sdB systems that possibly 
contain a neutron star or black hole companion: HE 0532-4503, 
HE 0929-0424, TONS 183 and PG 1232-136. 

These candidate sdB - neutron star (sdB-NS) systems were 
identified th rough spectroscopic observations of sdB binaries 
dGeier et al.L l2008i) . In sdB-NS systems, the optical spectrum 
is single lined and the compact companion cannot be detected 
directly optically. Assuming tidal synchronization between the 
orbital period and the sdB rotational period, it is possible to con- 
strain the binary inclination from measurements of the sdB sur- 
face gravity, projected rotational velocity and sdB mass. In their 
1201 Oi follow-up paper Geier et al. give a detailed description of 
this method. In binaries with periods shorter than 1 .2 days the 
orbit and stellar rotation are found to be synchronized. Thus, for 
such systems the companion mass can be derived from the sdB 
mass, whic h is known either fro m independent mass measure- 
ments (e.g. 'Fontaine et al.', '2008b or infe rred from binary pop- 
ulation synthesis (Han et al., 20^ l2003h . Our sources all have 
orbital periods shorter than these 1.2 days. In Table [T] we de- 
scribe further target parameters such as galactic coordinates, lo- 
cal background sky temperature, distance and estimated disper- 
sion measure. 

3. Observations and data reduction 

On 2008 Oct 15, 18 and 20 we observed these four sdB bina- 
ries with the Robert C. Byrd Green Bank Telescope (GBT). We 
used the lowest-frequency receiver PFl, which provided a band- 
width A/ = 50 MH z centered around 35 MHz. Every 81.92//S 
the Spigot backend ('K aplan et al.L l2005h recorded 2048 spectral 
channels in 16-bit total intensity. Over the three sessions, three 
test pulsars (Table |2]i and four sdB-star binaries (Table |3]l were 
observed. Due to scheduling and weather constraints, integra- 
tion times of the binaries varied, ranging from 12 to 68 minutes 
(Tables. 



Object 


hnt 


P 


DM 


Peak 5 /A' 




(s) 


(ms) 


(pc cm"^) 




PSRJ0034-0534 


60 


1.877 


13.8 


22 


PSRB0450-18 


30 


549.0 


39.9 


115 


PSRB1257+12 


30 


6.219 


10.2 


15 



Table 2. List of the previously known pulsars that were used as a 
check of the telescope, back-end and pulsar-search systems. Test 
source PSR J0034-053 4 is in a short-orb it 1.6-day binary with a 
0.16 M0 companion (IBailes etami994 . 



For each system in our sample, we derived the expected 
dispersion measures (DM) from their measured distance, us - 
ing the NE 2001 free electron model dCordes & Laziol |2002|) . 
These expected DMs, listed in Table [1] are all well below 
lOOpccm"^. Up to 200pccm"^ the trial DMs were spaced 
such that our search remains sensitive to pulsars with rota- 
tional peri ods down to 1 m s ; faste r pulsars are not generally 
expected dChakrabartv et al.l l2003h . Beyond 200pccm"^ the 
intra-channel smearing begins to dominate, so trial DMs up to 
lOOOpc cm""* where created at a lower time resolution. We next 
searched these trial DMs for pulsar signals, in both the time and 
frequen cy domain, using the PRESTO data reduction package 
(iRansomlBooil) . In the time domain we searched for single dis- 
persed pulses of radio emission. In the frequency-domain search 
we have to take into account the acceleration present in these 
binary systems. One could potentially derive the neutron-star or- 
bital phase from the optical modulation of the sdB star, and thus 
estimate the acceleration at the time of the observations. This 
however relies on the assumption that there is no lag between 
the optical light curve and the orbital motion. We have taken the 
conservative approach to search over the full range of possible 
accelerations throughout the orbit. As each integration time f,„, 
was shorter than 10% of the binary peri od, our search in pe- 
riod and period-derivative space sufficed (IJohnston & Kulkarnil 
1991). Using standard PRESTO routines, for each observation 
the radio-frequency interference was flagged, candidate period 
signals were sifted to remove harmonics, and the top 30 can- 
didates we folded and further refined by searching nearby dis- 
persion measures, periods and period derivatives to maximize 
signal-to-noise ratio. All candidate plots were subsequently vi- 
sually inspected. 

At the start of each observing session a known pulsar was 
observed; as listed in Table |2] these were all blindly re-detected 
by the above sdB-NS search pipeline, thus confirming the efifec- 
tiveness of the search method. 

4. Results 

No new pulsars were detected toward any of the four sdB stars 
in our sample. 

We next investigate how strongly these non-detections rule 
out the presence of an MSP. We derive the minimum detectable 
mean flux density S min for each observation from the pulsa r ra- 
diometer equation dDewev et al.L 119851; iBhattacharvai 1 1 9981) : 

_ {SIN,„i„)T,y, j W 

In this equation S IN,„m is the minimum signal-to-noise ratio 
at which a pulse profile can clearly be recognized as a pulsar pro- 
file. The system temperature r,,., is defined as T^ys - Tsys,GBT + 
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Object 


I 


b 




^orb 




d 


DM 








(Mq) 


(d) 


(K) 


(kpc) 


(pc cm"'') 


HE 0532-4503 


251.01 


-32.13 


1.4-3.6 


0.2656 + 0.0001 


17.0 


2.8 


44 


HE 0929-0424 


238.52 


32.35 


0.6 - 2.4 


0.4400 + 0.0002 


15.0 


1.9 


38 


PG 1232-136 


296.98 


48.76 


2.0 - 7.0 


0.3630 + 0.0003 


21.9 


0.57 


14 


TONS 183 


287.22 


-83.12 


0.6 - 2.4 


0.8277 + 0.0002 


19.7 


0.54 


12 



Table 1. Pulsar se arch targets based on sdB-NS candidate binaries identified bv lGeier et alJ (l2008h . Compa nion mass and binar y 
period taken from iGeier et al.l (l2008|) and references th ereii i. Sky tempera t ure at 408 MHz extracted from iHaslam et al] (Il982h . 
Distance determinations taken from Liske ret al.l (l2005h and lAltmann et al.) ( |2004|) . Dispersion measure estimate based on the NE 
2001 electron model dCordes & Lazia, 20021) . 



Tsky where Tsvs.gbt - 46 k Q. For each source, Ts^y was extracted 
from lHaslam et al.1 (Il982h and scaled from 408 to 350 MHz us- 
ing the Tsky °^ y'^'^ scaling law from Lawson et al. ( 1987). The 
gain G for the GBT is 2 KJy"' (iPrestage etalil2006h . the num- 
ber of polarizations Hp, added for these total intensity observa- 
tions, is 2. Integration time f,„, and bandwidth A/ for each obser- 
vation are described in Section [3] Finally, W is the pulse width 
and P is the rotational period. As these are unknown for non- 
detections, we use the average pulse duty cycle W/P=0.12+0.10 
determined by averaging over the MSP Wsn/P du ty cycles in the 
ATNF pulsar database (Manches ter et al.L l2005h . Here and be- 
low, we define MSPs in the ATNF database as those sources with 
periods shorter than 50 ms (thus, moderately to fully recycled) 
and magnetic fields lower than 10" Gauss. For each source, we 
list the resulting S,„i„ value in Table [3] 

The errors on our values of S,„m have contributions from the 
error on the average pulse duty cycle and from the systematic 
eiTor on the minimum signal-to-noise ratio S IN,nm- The S /N,„i„ 
at which a pulsar can be detected depends on the shape of its 
profile: a strongly peaked profile at high DM is more easily iden- 
tified as a pulsar than a same-S/N sinusoidal profile at low DM, 
as the latter can also be terrestrial interference. We estimate this 
systematic error as follows. To set a lower limit to the possi- 
ble iS INmin) range we inspected the pulse profiles of the test 
pulsars for increasing fractions of the total integration time. We 
concluded that starting from a peak {S INmin) of 7 the pulsar re- 
detections become unambiguous. As an upper limit to the range 
of {S /N,„in) values we adopt a value of 10. 

To estimate the completeness of our search we derived 
pseudo luminosity L-S cfi limits for each of the binaries in our 
sample and compared them with the pseudo luminosities of 
known MSPs. We define the completeness of our pulsar search 
as the percentage of known millisecond pulsars that have pseudo 
luminosities higher than our derived pseudo luminosity upper 
limits, i.e., as the percentage of known MSPs that would be de- 
tected if placed at the distance of the candidate. We searched 
the literature for distances d to the binaries in our sample and 
used the most recent values. For HE 0532-4503 and HE 0929- 
0424 these are from Lisker et al. (2005), while the dist a nces t o 
TONS 183 and PG 1232-136 are from lAltmann et all (|2004 . 
where TON S 183 is known as SB 410. The former article claims 
an error on the distance of 10%. In the latter article distance er- 
rors are not estimated, so for those sources we also propagate a 
10% error All distances were derived from models of sdB atmo- 
spheres combined with the magnitude measurement. 

In Figure[T]we compare our limits to the pseudo luminosities 
of the millisecond pulsar population. For this statistical compar- 
ison, we have selected all 5 MSPs with k nown 400 MHz fluxes 
in the ATNF database ( Manchester et al.L 12005 ). We scale these 



Object 






V . 

^ nun 


L 


C 




(s) 




(mJy) 


(mJy kpc^) 


(%) 


HE 0532-4503 


720 


0.41 


±0.21 


3.2 ± 1.8 


88 


HE 0929-0424 


900 


0.35 


±0.18 


1.3 ±0.7 


96 


PG 1232-136 


1080 


0.37 


±0.19 


0.12 ±0.07 


100 


TONS 183 


4080 


0.18 


±0.09 


0.05 ± 0.03 


100 



Table 3. Integration time f„„, flux upper limit S,„m derived from 
radiometer equation, pseudo luminosity limit L calculated from 
the flux upper limit, and completeness C compared to the known 
MSPs in the ATNF database. 




http : // sci ence . nrao . edu/gbt/ obsprop/GBTpg . pd£ 



12 3 

Pseudoluminosity ( log(mJy kpc2) ) 

Fig. 1. Cumulative histogram of pseudo luminosities for the 
known MSP population in the ATNF database. The left axis is 
labelled with the cumulative fraction of MSPs; the right axis is 
labelled with survey completeness for each candidate. The er- 
rors on the pseudo luminosity upper limits contain eiTors on the 
distance, pulse duty cycle and our estimate of the relevant range 
of (S INmin) values used with the pulsar radiometer equation (see 
Sectionllli. 



400 MHz fluxes to our central observing fr equency of 35 MHz , 
using the -1.6 average spectral index that iKramer et al.l (Il998h 
find for MSPs. We find that our limits exclude with high con- 
fidence the presence of pulsar emission: in Table [3] we list the 
exact values for this completeness C. 

5. Discussion 

As no radio pulsations were found in any of our observations, we 
place upper limits on the pulsed radio emission from the putative 
neutron stars in these systems. Assuming the known MSP lumi- 
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nosity distribution, we can exclude the presence of such pulsed 
radio emission with 100% certainty for 2 systems: PG 1232-136 
and TONS 183 (Figure[B. For 2 other systems, HE 0532-4503 
and HE 0929-0424, we were not able to put as strict upper lim- 
its on the pseudo luminosity due to their larger distance and the 
shorter integration time for these pointings. For HE 0929-0424, 
our survey would have detected 96% of MSPs at the HE 0929- 
0424 distance. For HE 0532-4503, this completeness C is 88% 
(Table|3] Figure©. 

In relation to these non-detections we now discuss two se- 
lection effects: the fraction of the sdB lifetime during which an 
MSP is on; and the beaming fraction, which is the fraction of sky 
over which an MSP beam sweeps, and thus from which the MSP 
is in principle detectable. As outlined in SectionlT] the MSP and 
sdB star are formed simultaneously; as MSPs have ages up of to 
10^"'" years and sdBs have ages up to 10^ years, an MSP formed 
in a binary with an sdB star will shine for the entire age of the 
sdB, and longer. We thus conclude there is no age bias against 
detecting MSPs around sdB stars. There is however a non-zero 
chance that for some of our non-detections a bright MSP emis- 
sion beam is present, bu t missing the Earth. The beaming frac- 
tion of MSPs is 0.7±0.2 jKramer et al.Lll99 8'); high compared to 
normal pulsars, but not unity and thus an important factor. 

So far we have discussed MSPs, as any regular (non- 
recycled) pulsar originally present in these systems will have 
long shut off before the formation of the sdB began (cf. Section 
[T]l- For any non-recycled pulsars that are in the beam by chance, 
the S min in Table [3] list the minimum detectable flux for that 
pointing. The pseudo-luminosity distributions of non-recycled 
and millisecond pulsars are similar ( Manchester et al.„,2005) . but 
as the distance to any chance-coincident non-recycled pulsar is 
unknown, their pseudo luminosity is unknown too. 

We now com pare the res u lts of our search, triggered by 
the sample from iGeier et al ] (|2008), with the follow-up re- 
sults of furt her spectroscopic investigation recently published in 
iGeier et al ](^010). For PG 1232-136, where we put a very strict 
(100% complete) limit on the pseudo luminosity, iGeier et aH 
( 20101) now report the mass of the unseen companion Mcomp 
to be higher than 6 Mq. This indicates the compact compan- 
ion is a black hole; therefore no pulsed radio emission is ex- 
pected from that system, in agreement with our non-detection. 
For TON S 183 the new spectroscopic results point to a Mcmp 
of 0.9 Mq but the error bars allow for a low mass white dwarf 
(WD). Our non-detection of pulsed radio emission is compat- 
ible with a WD companion, and since the constraint we put 
on pulsed radio emission from TON S 183 is strict (100% com- 
plete) we rule out an MSP beamed towards Earth. For HE 0929- 
0424 a Mcomp slightly above the Chandrasekhar limit is reported, 
but with error bars that also allow for a high mass WD. Our 
96%-complete non-detection in this case cannot exclude either 
possibility but an MSP is unlikely. The system with the least 
strict pseudo luminosity upper limit, HE 0532-4503, is reported 
to contain a companion to the sdB star of 3 Mq. Even if one 
allows for an sdB star as light as 0.3 Mq the companion re- 
mains more massive than the Chandrasekhar limit. Given our 
non-detection corresponding to the 12th percentile of the known 
MSP population (88% complete), an MSP is not likely, but this 
system remains an interesting candidate for deeper radio follow- 
up. Further observati ons of this source and of new candidate 
sdB-NS systems from lGeier et alJ (1201 Ol) are currently ongoing. 



6. Conclusions 

We have searched for pulsed emission from potential pulsar 
companions of four subdwarf B stars. No pulsed emission was 
found down to luminosities corresponding to, on average over 
the 4 sources, the 4th percentile of the known millisecond pulsar 
population. 
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